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Abstract 

Recently a local true (completely gauge fixed) Hamiltonian for spherically symmetric collapse was 
derived in terms of Ashtekar variables. We show that such a local Hamiltonian follows directly from 
the geometrodynamics of gravity theories that obey a Birkhoff theorem and possess a mass function 
that is constant on the constraint surface in vacuum. In addition to clarifying the geometrical 
content, our approach has the advantage that it can be directly applied to a large class of spherically 
symmetric and 2D gravity theories, including p-th. order Lovelock gravity in D dimensions. The 
resulting expression for the true local Hamiltonian is universal and remarkably simple in form. 

PACS numbers: 04.60.Ds 
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I. INTRODUCTION 



Despite considerable progress in recent years in string theory and loop quantum 
gravityj2], relatively little is known about the quantum dynamics of self-gravitating systems. 
General covariance provides the classical theory with much of its underlying elegance but it 
also leads to deep conceptual problems due to the need to fix a gauge (i.e. choose a coordi- 
nate system) and the associated problem of timej^. The situation is further complicated by 
predictions, primarily from string theory, of quantum induced, higher order curvature terms 
in the action and the possible existence of more than three spatial dimensions. 

Since the groundbreaking work of Berger et. a/.j^ and Unruhjs] in the 1970's spherically 
symmetric scalar fields have provided a fruitful testing ground for the classical and quantum 
dynamics of gravitational collapse. For example, as first shown by Choptuikjol, spherical 
collapse exhibits interesting critical behaviour. The Birkhoff theorem, which states in essence 
that there is no spherically symmetric gravitational radiation, is a key simplifying feature, 
since there are at most only a finite number of gravitational dynamical modes to consider. 
At the quantum level, the Birkhoff theorem guarantees that the vacuum theory can be 
quantized exactly, and this has been done via many different techniques j?]. 

Despite the drastic simplification afforded by the Birkhoff theorem, the dynamics of the 
matter degrees of freedom is greatly complicated by gravitational self-interactions. The 
study of the quantum dynamics has to date been rendered intractable by the fact that 
completely fixing the gauge appears to invariably lead to a non-local Hamiltonian for the 
1 15[. (For a current discussion see |8|, l9|.) 



10|, who derived a local 



matter field 

A significant step has recently been made by Alvarez et. al. 
gauge fixed Hamiltonian to describe the spherically symmetric collapse of a massless scalar 
field coupled to Einstein gravity in four spacetime dimensions. They described the gravita- 
tional field via Ashtekar's connection variables as presented most recently by Bojowald and 
Swiderski H ] . 

In the following, we prove that geometrodynamical arguments can be used to obtain a 
similar, but algebraically simpler, local Hamiltonian in a large class of gravity theories that 
obey a Birkhoff theorem. The suggestion that the analysis is greatly simplified in terms of 
the standard ADM variables was made independently by V. Husain[12] and W. Unruhjisl 
and verified explicitly by Unruh for 4D Einstein gravity. We are able to achieve a further 
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simplication by parameterizing the gravitational variables using the mass function and its 
conjugate. In addition to allowing us to consider simultaneously a large class of theories, 
our choice of parameterization and gauge yields a reduced Hamiltonian that is very simple 
in form. This in turn facilitates the discussion of the consistency of the gauge condition for 
generic initial data. 

The theories to which our analysis applies include generic 2D dilaton gravity[l^ and p-th 



order Lovelock gravity 



15, 



16| in D dimensions |29j. Lovelock gravity theories were originally 



constructed as the most general manifestly covariant higher dimensional versions of Einstein 
gravity for which the equations of motion are second order in the time coordinate- this 
ensures reasonable causal behaviour for the propagation of gravitational fields [15]. There 
are string theoretic arguments that point to Lovelock gravity as the low energy limit of 
consistency conditions for the propogation of strings in curved spacetimes [l^]. In particular, 
these stringy consistency conditions contain higher order curvature terms in the same specific 
linear combinations that occur in Lovelock gravity. If higher dimensional black holes can be 
probed via the LHC, then it may be possible to rule out certain terms in Lovelock gravity 
theories 18|. Finally, and most importantly for what follows. Lovelock theories are known 
to obey a generalized Birkhoff theorem fl9. 120I. 

The paper is organized as follows: in the next Section we describe the geometrodynamics 
and Hamiltonian framework of theories of gravity that obey a Birkhoff theorem and hence 
possess a mass function^]. Section III shows in full generality how to gauge fix in such 
theories so as to obtain a true local Hamiltonian. Section IV applies the general framework 
to p-th order Lovelock gravity in D dimensions and derives the local Hamiltonian while 
the following Section discusses boundary conditions. Finally, we close with a summary and 
prospects for future work. 



II. GEOMETRODYNAMICS IN GENERAL 



We will work with the following ADM decomposition using Kuchar's 22|] parameters and 
conventions: 

ds^ = -N^dt^ + A^dr + Nrdtf + R^dVl''-^ (1) 

The metric functions A^, A, Nr, R are arbitrary functions of the two spatial coordinates (r, t), 
while R{r, t) is the areal radius of the spherically invariant {D — 2) spheres. The following 
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analysis also applies directly to gravity in two space-time dimensions where the role of R 
is played by a scalar field (the dilaton) non-minimally coupled to the curvature. Standard 
Hamiltonian techniques reveal that for the spherically symmetric ansatz above any coordi- 
nate invariant gravitational Lagrangian can be written, up to boundary terms, in the form: 



L(^) = PaA + PrR - NH^^^ - NrHf\ (2) 

where H^^'^ and hP are the gravitational parts of the Hamiltonian and diffeomorphism 
constraints, respectively. The transformation properties of the ADM variables under changes 
of spatial coordinates require the gravitational diffeomorphism constraint to be: 

iff) = -P^A + PrB!. (3) 

Adding a minimally coupled scalar field ip yields the Lagrangian: 

Ltot = L^""^ + L^"'^ = jdr (^P^A + PrR + U^ij - Ar(i/(^) + ^) - Ar,(i/f ) + , 

(4) 

where 

P^'"^ - Ul^^ + W V (5) 



2 \R^-^ 

is the energy density of the scalar field without self-interactions. This could be generalized 
to other types of matter subject to the restriction that it does not contain derivatives of the 
metric. The form of the matter part of the diffeomorphism constraint is again dictated by 
the transformation properties of the fields. For a minimally coupled scalar field iJj: 

Hi^^ = li^ij'. (6) 



In what follows we show that a local Hamiltonian as in 10| can be derived using Kuchar 



geometro dynamics 24] in a large class of theories that satisfy a Birkhoff theorem, i.e. for 
which the spherically symmetric vacuum solution is static and is parametrized by a mass 
function M that is constant on the constraint surface. In a fully dynamical setting any 
spherically symmetric metric can be written in Schwarzschild-type coordinates as: 

dr^ = -F{R; M)N^dT^ + F-\R; M)dR^ + R^dQ^~^ (7) 



where M = M{R, T) is the (generalized) Misner-Sharp mass function 2l| and Ns = Ns{R, T) 
is the lapse function. 



In D dimensional spherically symmetric Einstein gravity the Misner-Sharp mass can be 



written 



2l|: 



M 



R 



D-3 



:i - \VR\ 



R 



D-3 



whereas in general Lovelock gravity it is of the form 16| : 

M := > ' a,„)R''-^-^f{l - FY 



[n/2] 

p 



(8) 



(9) 



The mass function M and the areal radius R are invariant under coordinate transformations 
preserving the spherically symmetric form of the metric. 

The essence of geometrodynamics||22i] is the construction of the canonical transfor- 
mation relating the ADM variables {A, Pa), {R, Pr) to the new phase space variables: 
{Pm,MY{R,Pr)- This enables one to extract the physical/geometrical properties of the 



theory 
al 



23|. The geometrodynamics of 4D Gauss-Bonnet was first described by Louko et. 



24| and has recently been completed for p-th order Lovelock gravity 25|, |26|]. Previous 



analyses considered only the static, vacuum case, but it also can be shown to go through for 
the dynamical metric ([7]) above [26|. 
It is useful to define 

y:=N'\R-NrR!), (10) 

n 

since R can only appear in the action in this invariant combination [22]. Doing the coordinate 
transformation from {R,T) to (r, t) in (JTj) and matching coefficients of dr'^, etc. in ([T]) one 
obtains the key relations: 



NsT' 

F{R-M) = j^-r 



Ay 
F ' 



As shown originally by Kuchar [22'] for vacuum spherically symmetric gravity: 

Pm := = -N.T\ 



Ay 
' F 



(11) 
(12) 

(13) 



is canonically conjugate to the mass function M. The analysis is identical in the dynamical 
case, since appears only in the relationship between Pm and the Schwarzschild time T. 



Kuchar 



221 ] argued, based on the form of the diffeomorphism constraints before and after 



the transformation, that Pr will obey: 

Pr = Pr- P'aA/R! - PmM'/R'. 



(14) 



After the change of variables from (A, Pa, -R, Pr) to (M, Pm, R, Pr), Eq.(|2]) can be written: 

L^G) ^ p^^^j ^ p^p _ p^ff(G) _ p^^Hp. (15) 

A sufficient condition for Kuchar' geometrodynamics to work is that the gravitational 
action depend linearly on A. This is the case for generic 2D dilaton gravity and for generic 
Lovelock theory [26]. Following [25'] we write the gravitational lagrangian as : 

L(^) = B,{R, y, A) A + Bo{R, y, A). (16) 

Clearly 

PA = B^iR,y,A), (17) 
can be inverted to solve for y purely in terms of A, Pa, R- That is 

y = y{A,P^,R). (18) 

The crucial property of the above is that y does not depend on Pr so that we can use ( fT3|l 
and f|T2|) to solve algebraically for A and Pa in terms of M, Pm and R. The Hamiltonian 
derivable from (fT6|l then takes the form: 

n^^^ = PaA + PrR-L^'''' =PRR-Bo{R,y,A) 

= NrPnR' + NPnyiA, Pa, R) - Po(A, Pa, P), (19) 

where we have used the definition of y to get the last line. Thus the Hamiltonian is linear 
in Pr. The first term will of course contribute to the diffeomorphism constraint while the 
second will appear in the Hamiltonian constraint. This will be important in what follows. 

We next assume the existence of a Birkhoff theorem, so that in vacuum the solution can 
be put in Schwarzschild form ([71) with M = constant. In this case, there must be a linear 
combination of the gravitational Hamiltonian and diffeomorphism constraints such that: 

^if (G) + 5^^(G) ^ ^20) 

so that on the constraint surface M' = 0. The coefficients depend on the phase space 
parameters but must be non-degenerate almost everywhere on the constraint surface. By 
virtue of (fT2l) and (fT9|) . respectively, M' and the coefficients A,B are independent of Pr. 
They do, however, depend on (A , Pa) and hence on (M , Pm)- 
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Using (l20l) the full Hamiltonian in the presence of matter can be written: 
Htot = ^ ( - BHf^ + j + Nr{PRR' + PmM' + H^V^') 

= ^ f -M' + 5n^^' + + f iV, - j {PrR! + PmM' + H^V^')- (21) 

One finally gets the following very general and intuitive result for the Hamiltonian constraint: 

H = -M' + ^p^"^) + Bli^^' ^ 0. (22) 

It is intuitive because it says that on the constraint surface the gradient of the mass function 
equals the energy density of the matter fields (if^"*)) plus a gravitational self interaction term 
9]. The above is a more general version of the result in jol for generic 2D dilaton gravity, 
which includes Einstein gravity with cosmological constant in any dimension. 

III. GAUGE FIXING AND THE LOCAL HAMILTONIAN 

Without loss of generality, choose the spatial coordinate to depend only on the areal 
radius, so that R = R{r). The consistency condition R = then determines the shift via 
the total Hamiltonian ( 12T]1 : 

Nr = N^. (23) 



The partially reduced Lagrangian is then: 



-'pr 



dr [ PmM + H^^ - ) , (24) 



where H in fl22|) depends only on {ip, H^, M, Pm)- 

Given the above generic geometrodynamic structure we can now complete the gauge 
fixing by choosing the mass function as follows: 

X = M-/(r,t)^0. (25) 

for some suitable function f{r,t). Note that a more general choice, of the form: 

Mir,t)=gii;)fir,t), (26) 



for arbitrary g{ip), as done in 10|, is also possible. (!26|) can be implemented by performing 



the following canonical transformation in 

M = M/gii;); Pj^, = g{i;)PM; 



ij = i^; H^ = H^ + PmM^, (27) 



and then making the gauge choice M = f{r,t). The only change is that H in becomes 
a more complicated function of P^^ which in turn makes the algebra messier. We henceforth 
restrict consideration to the simpler case (l25l) . 

It is important to note that ( l25l) does not provide a good gauge fixing condition in 
vacuum since M commutes weakly with the gravitational part of the Hamiltonian constraint. 
However when the matter fields are not zero, M does not have zero Poisson bracket with the 
Hamiltonian constraint by virtue of the latter's dependence on through the coefficients 
A and B. In this case, one can use the consistency condition x to determine the 
lapse and the Hamiltonian constraint H ^ determines Pm in terms of the matter phase 
space variables ip, H^ and the (gauge fixed) mass function M. The fully reduced Lagrangian 
density is therefore: 

Lred = H^,^ + Pm(^, H^)M, (28) 
with corresponding reduced Hamiltonian density: 

Hred= -PM{lp,Ili,)M (29) 

By inspection of H it is clear that the reduced Hamiltonian will be a local, algebraic, time 
dependent function of the matter fields. Here and in the following we use M(r, t) to denote 
the gauge fixed mass function, so that it is no longer a phase space variable. 

In order to make the above construction explicit and verify its broad range of applicability, 
we now turn to generic Lovelock gravity in D dimensions. 



IV. P-TH ORDER LOVELOCK GRAVITY IN D DIMENSIONS 



The action for generic p-th order Lovelock gravity in D dimensions can be written 15|: 

/ = ^ / rf"a:v/=^ 5^ (30) 

where S^lZpp '■= ^[p-^■■■^pp\ ■ Each is a coupling constant of dimension (length) ^{p-i)^ with 
0(0) proportional to the cosmological constant. The p = 1 term gives the Einstein- Hilbert 
contribution. The specific combination of curvature terms in each of the £(p) guarantees that 
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the equations for Lovelock gravity are second order in the metric components and ghos t-free. 



19 



20| 



As mentioned above the generic theory obeys a generahzed Birkhoff theorem 
For Lovelock gravity, it has been shown [26] that (12 Op is satisfied with 

where 

K' = j{R'^-F'Pl). (33) 

as implied by (fT2|) and ( IT3|) . After setting R = r, the consistency condition R = determines 
the shift to be: 

Nr = PmF-. (34) 
The partially reduced action is given in terms of geometrodynamic variables by: 

dr I PmM + - (-M'A^ + p(-) + PMFYi^i;') \ . 
For completeness we write down the equations of motion for the scalar field: 

A^ / 



(35) 



At this stage one can choose a variety of different gauges: Pm = yields Schwarzschild-type 
coordinates, while A = 1 yields fiat slice or Painleve-GuUstrand coordinates 27||. In order to 
obtain a local Hamiltonian we now specify the mass function as in ( j25l) . The consistency 
condition % for this gauge choice determines the lapse function: 

N M 



(37) 



A 2M'PmF + Fli^ij' 
The gauge condition is only valid if H} ^ 0, which, for A^ 7^ reads: 

2M'PmF + Fn^^' ^ 0. (38) 

As long as the above inequality holds, one can set the gauge fixing condition and the hamil- 
tonian constraint strongly to zero. The Dirac brackets in this case equal the original Poisson 
brackets since both ip and 11^ commute with the gauge fixing condition in this parametriza- 
tion. This yields the true Hamiltonian in (129|) where Pm is a local function of (^, 11^) 
obtained by solving the Hamiltonian constraint: 

M' 

— (1 - F'Pl) = p('") + PmFU^^'. (39) 



The solution to this quadratic equation for PmF is: 



— mf ■ 

For the classical theory, we consider only the Hamiltonian density with the positive root of 
the discriminant. Indeed, as we will see in the next Section this gives the result that in the 
case when p*^™-' — M'/F is small, one can define a perturbative regime in which, up to an 
additive constant, the Hamiltonian density ~ p*^™"-*, . 

The consistency condition f l38|) implies that discriminant A under the square root in pUj) 
not vanish. Given that we want the geometrical variables to be real, we require the stronger 
condition: 



< A := (n^^')' - ^(P^™^ - M'/F) 

r 



4 \RD-'^ 
The above is of the form 



4 V R^-^ 



(41) 



\A^ + -2M' /F\>\A^ - B\ (42) 

with = 27ji7b2 and B"^ = R^~'^ip'^/2, so that 11^^' = 2AB. Clearly a sufficient condition 
for this is that 2M'/F < 0. More generally, for a given set of initial data, at t = 0, say, 
it is relatively straightforward to implement a valid gauge condition by making a judicious 
choice of M'/F\t=Q = P{r)^0\. One possibility, which will be relevant for what follows, is 
simply: 

M' 

= p('")+e(r)/2 = /lY52 + e(r)/2, (43) 

^ t=o 

where e(r) is everywhere positive. The consistency then reduces to: 



\A^ + B^ + e\-\A^-B^\>0, (44) 

which is identically satisfied. We therefore claim that it is in principle possible to find a 
suitable gauge fixing condition of the form (123]) for generic initial data. The term "suitable" 
refers to the validity of the gauge choice fl55]) and the positivity of the square root at time 
t = 0. This does not address the issue of boundary conditions nor the possibility that the 
gauge may break down after some finite time t. Boundary conditions will be discussed in 
the next Section. 
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The resulting reduced action, given explicitly in terms of (V',!!^), is remarkably simple 
and universal in form: 



M M' 



dr H^ip - 



M' F 
M 



2M' 



AM' 



(45) 



This is the main result of our paper. 

We note that knowledge of the scalar field and its conjugate on a given slice allows the 
metric functions A^, Nr and in terms of the scalar field and its conjugate to be determined 
using fl371) . flMl) and fl33|) . respectively: Thus the complete spacetime can in principle be 
reconstructed from the data {ip, 11^) on each spatial slice, as long as none of the expressions 
diverge or become imaginary. Since, as argued above, the gauge choice fl2^ is invalid in 
vacuum, i.e. when the matter field and its conjugate vanish, we now turn to the important 
issue of boundary conditions. 



V. BOUNDARY CONDITIONS 



For the matter field, the usual boundary conditions consistent with asymptotic flatness 
in D space-time dimensions are usually chosen to be: 

ip — )• Br 2 2 j 

n^^Cr^-t, (46) 

so that 

^ CBr-\ (47) 

The only boundary terms that arise in the variation of fH5|) come from the variation of ip' . 
For e > 1, these will be finite as long as A 7^ 0, as required by the consistency condition on 
the gauge choice, and M /M' is regular at the boundary. 

Since the gauge fixing condition is problematic in vacuum, one might expect correspond- 
ing subtleties to arise in asymptotically fiat space times. These subtleties are revealed by 
inspection of of the metric functions. In fact (|39|1 implies that if p*^"^^ — )■ faster than M', 
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then FPm must go to unity. This in turn yields — )■ 0, which signals a breakdown in the 
coordinates. In addition M' — > faster than the matter fields is ruled out by f l39|) since it 
would required |1 — -F^P|j| > 1, which would result in < 0. The question then arises: for 
what physically relevant calculations is the local Hamiltonian useful? 

One possibility is to avoid asymptotically fiat spacetimes and explore more general sit- 
uations, such as those containing ingoing/outgoing radiation. Alternatively, one can note 
that the perturbative regime (for which the scalar obeys the usual free field equations) 
corresponds to: 

pM _ M'/F -> 0, (48) 
as follows from f l45|) . In this case, for the + sign in fH5l) . we find 

MM' , Mp^-M'/F , , 

Since Pu this is consistent with the expression for N j K above, i.e. in this limit 

N M 
A - 

Moreover, in this limit iV^ — ?• and A^ — )■ So if we can choose 

y ^ n^V^', (51) 

everything appears consistent. 

This suggests the following procedure for fixing initial data in calculations of spherical 
self-gravitating collapse: First assume that the matter obeys the fall off conditions at 
t = 0. Then choose gauge fixing to be of the form: 

M = a(r) + b{r)t. (52) 

Finally determine the functions a(r) and b{r) so that on the initial slice: 

n^^' + 0(r-"). (53) 

Note that this choice of M(r, 0) matches (143|) in the previous section. One can then evolve 
the matter field according to the equations of motion ( 136|1 until the gauge fixing condition 
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M' 






a'(r) 


~y 


t=0 


1 


— a{r)/r 


M 






b{r) 


T 


t=o 


1 


— a{r)/r 



breaks down or one of the metric functions diverge. Note that from the equations of motion 
for the matter field, the asymptotic time evolution for the matter field is governed by: 

This suggest thats the asympotic relationship between matter and gauge fixing will be 
preserved for a finite time. 



VI. CONCLUSION 



We have shown that a local Hamiltonian can be obtained for spherically symmetric mass- 
less scalar field dynamics in any theory of gravity to which the geometrodynamics of Kuchar 
can be applied. The class of theories for which this holds includes the generic Lovelock case 
as well as generic 2D dilaton gravity. The reduced action ( H5l) is remarkably simple and 
universal for all Lovelock theories in any dimension. The existence of a local Hamiltonian is 
a direct consequence of Eq. f|T9|) which shows that the gravitational Hamiltonian is linear in 
Pr. In fact there are two terms, one each of H'^'^^ and Hi'^\ The coefficients A and B that 
determine f l39|) are, up to an overall normalization factor, precisely those required to cancel 
the two linear terms: AR' + By = 0. The overall normalization factor should be the same 
for all such theories, and is such as to yield NA/R' as the coefficient of M' in fl^T]) . This 
coefficient plays a role in determining the boundary conditions, since the boundary term 
in the Hamiltonian comes from the variation of M'. We conjecture that the form of the 
true Hamiltonian f HOj) in this type of gauge is the same for all theories that obey a Birkhoff 
theorem and possess a mass function M that is constant on the constraint surface in vac- 
uum. The gravitational back reaction is determined completely by the dependence of F on 
M, as one may expect from physical grounds. The extension to theories with cosmological 
constant is straightforward, since it merely modifies the relationship between F and M. 

We note that our geometrodynamic formulation and family of gauge choices results in a 
reduced Hamiltonian that is much simpler than that of [10]. This makes it easier to discuss 
both the advantages and potential disadvantages of this family of gauges. Specifically, we 
were able to show in Section IV that the reality condition on the square root is directly tied 
to the validity of the gauge choice and more importantly, we were able to find consistent 
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gauge choices for generic initial data. 

Despite the universahty and elegance of the final reduced Hamiltonian, there are potential 
problems in implementing this gauge fixing procedure in asymptotically flat spacetimes, 
related to the fact that the gauge condition is not valid in vacuum. We have proposed 
one possible approach to this problem for the investigation of spherical collapse, but the 
viability of this approach requires further investigation. As mentioned in Section III, another 
possibility is to consider more complicated gauge fixing conditions that depend on both ijj 
and M, as done in [10]. 

Finally, we note another potential application for the local Hamiltonian, namely the 
calculation of scattering amplitudes for spherically symmetric self-gravitating matter using 
the phase space path integral. In principle, the relatively simple algebraic form of the local 
Hamiltonian might provide simplifications over the method used in {28 1, for example, to 
calculate the tree level four point interaction. However, it should be noted that in the case 
of Hamiltonians quadratic in the momenta it is possible to perform the integral over the 
momentum to obtain a simple covariant path integral with a Lagrangian that is quadratic 
in the velocities. In the present case, the integration over the momentum is not trivial 
and given the non-polynomial nature of the Hamiltonian, it will yield an expression that 
is a non-polynomial function of the velocity. It is nonetheless an intriguing question that 
deserves further attention. 
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